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Abstract 



^NJ ■ The decays Kl — » 77 and i^x, — ► ^ + £~ 7 are studied at the leading order p 6 

> 

in Chiral Perturbation Theory. One-loop contributions stemming from the 

pg ■ odd intrinsic parity | AS |= 1 effective Lagrangian of order p 4 are included 

o ' 

C~-- ■ and shown to be of possible relevance. They affect the decay Kl — > 77 

(~| ■ adding to the usual pole terms a piece free of counterterm uncertainties. In 

a: 

q; the case of the Kl — » £ + £ 7 decays the dependence of the form factor on 

the dilepton invariant mass requires a counterterm. The form factor may re- 
V"! ■ ceive a sizeable contribution from chiral logarithms. Including considerations 

from the Kl — > it + 7t~ 7 direct emission amplitude, we obtain two consistent 
scenarios. In one scenario the long distance contributions from the one-loop 
terms are important, while in the other they are marginal. In both cases the 
counterterm is shown to be significant. 
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The radiative decays of the long lived neutral K meson, Kl, have been the subject 
of numerous theoretical analyses, in Chiral Perturbation Theory (xPT) |I[. Among these 
decays, the decay Kl —> 77 and the Dalitz decays Kl — > ^ + ^~7 still show open issues 
which we address in this paper. Both types of decays proceed at order p 6 in ;\TT. The 
pole-type terms given by 7r° and the rj 8 intermediate states |2|-f| are individually of leading 
order p 4 , but their added contribution cancels at this order due to the Gell-Mann-Okubo 
relation || and is in effect of order p 6 . It is, therefore, important to perform a complete 
analysis of the decays at this order. We show that the one-loop terms of order p 6 stemming 
from the odd-intrinsic parity terms of the | AS |= 1 non-leptonic effective Hamiltonian of 
order p 4 and not previously evaluated, affect both types of decays in a potentially significant 
way. We find, in particular, a contribution to the amplitude of the two-gamma decay 
proportional to M\ + M% that can be a substantial addition to the pole-type term. In 
the case of the Dalitz decays, the experimentally observed non-trivial form factor in terms 
of the dilepton invariant mass is usually described using the VMD model of Bergstrom, 
Masso and Singer ||, which gives a satisfactory parameterization. The work of Ko |5| 
provides a thorough study of VMD in radiative decays. Recently, other models like the 
factorization model (FM) and the weak deformation model (WDM) have been considered 
as well @. Our calculation includes long distance pieces given by the chiral loops, and a 
short distance piece due to a counterterm needed to renormalize the loops is determined 
by fitting to the data. For our purpose, the experimental situation is good. The two- 
gamma mode is experimentally very precisely known: BK(Kl — > 77)= (5.92 ± 0.15) x 10~ 4 
||. The branching ratios of the Dalitz decays are on the other hand known within 10%: 
BR(K L -> e+e"7) = (9.2 ± 0.5 ± 0.5) x 10~ 6 (NA31 and BNL-E845 collaborations), 
and BR(K L -> ll+ll-^) = (3.23 ± 0.23 ± 0.19) x 10~ 7 (E799 collaboration JTJ). The rate 
of the fi + fi~ mode is about 30% larger than it would be for the case of a point-like form 
factor. The form factor also shows up clearly in dilepton invariant mass distribution of both 
modes. Using this data and some input from the Kl — > 7r + 7r~7 decay [|T2|] , we are able to 
two phenomenologically acceptable scenarios. 



We work in the standard framework of %PT. The octet of pseudoscalar mesons is pa- 
rameterized by the SU(3) matrix 
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where Fq ~ F„- = 93 MeV is the pion decay constant in the chiral limit. For the processes 
we consider the only gauge couplings we need are those to the EM field. We thus denote: 



Ln-BitfDpU 
D lt U = d lt U-iA„[Q, U], 



(2) 



where Q is the quark charge matrix. 

We need the strong interaction effective chiral Lagrangian up to order p A |Tj| , including, of 
course, EM couplings and the Wess-Zumino (WZ) Lagrangian. The non-leptonic | AS* = 1 | 



effective chiral Lagrangian is also required up to order p 4 [14|. It contains an 8 and a 27 of 
SU(3); we advocate the AI =1/2 rule to disregard the marginal effects of the 27 piece. The 
octet piece of order p 2 is 



jC\ 2 2 sl=1 = ^-C 8 Tr(\ 6 L,IS), 



(3) 



where Cs = 3.12 x 10 is given by 

C 8 
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Here Sj, q are the sine and the cosine of CKM angles in the Kobayashi-Maskawa represen- 
tation, and g§ ~ 0.5 is an effective low energy coupling. At order p 4 there is a long list of 



terms (forty eight octet terms when external sources of different kinds are included) |H 
The relevant terms to us are those that affect the virtual transitions Kl — > 7r°, 7/g, includ- 
ing terms that break SU(3) symmetry due to the quark masses, and terms that give the 



amplitudes Kl — ► 7r + 7r 7 and Kl — > 7r + 7r 77. The first are contained in the terms bilinear 
in the meson fields of the 0(p 4 ) Lagrangian |H[ , and the second are given by odd-intrinsic 



parity terms of which only two combinations of the operators Of in (14| appear, namely, 
W29 = 0| 8 + 0| 9 , and W31 = Of + O^. The piece of the Lagrangian £| AS=1 i containing 



these two terms can be written as follows W2 



4§3i = ^i {a 2 Tr(A[U^F^U, L,L V )) + a 4 Tr(AL„) Tr{{FT ~ U^F^U)L V )} H-h.c. 

(5) 

where A = |(Aq— 1X7)- In our case, where the only external gauge field of relevance is the EM 
field, we have Fg = F£" = F^Q, and F^ = ^ upa F pu is its dual. Since the Lagrangians of 
order p 2 are of even intrinsic parity, they alone cannot generate odd intrinsic parity terms. 
Therefore, the low energy couplings a 2 and 04 are renormalization-scale independent. 

In the absence of CP violation Kl coincides with the CP-odd state K 2 = 4^(K° + K°), 
and the pieces of C 29 31 leading to the transitions K L — >■ 7r + 7r~7 and K L — ► 7r + 7r~77 are 
explicitly: 

87r 2 F(f 
+ a 4 dfj,K L (n~d u n + - n + d u n~)] d p A a , 
e 2 

- 2 a 4 d p K L A u }d p A a 7i + n- + ... (6) 

In the second Lagrangian we do not display terms that do not contribute in our calculation. 
The SU(3) singlet meson rji, being a massive state due to the U^(l) anomaly, can be 
integrated out or, equivalently, included explicitly in a meson pole dominance model || f| . 
The T]x contribution to the decays is naturally of order p 6 . Since 771 and 7/g mix with each 
other due to SU(3) breaking, we find it more convenient to include 771 explicitly. We need 
to consider both, the terms that contribute to the virtual transition Kl — > 771 as well as the 
WZ term, which leads to the standard amplitude for rji — ► 77. The mixing angle between 



the 7] mesons and the decay constants as obtained from their two-gamma decays and the 
C(p 4 ) result for F m [|IJ are: 

9= -21±5deg 

" r/8 1 ■ O .r jj- 

^ = 1.1^ (7) 

The virtual weak amplitude Kl — > rjx is not known, and is characterized in the following by 
a parameter k. 

In the limit of CP conservation, the Kl — > 77* amplitude has the most general form: 

A(K L - 77*) = F(t) iT" e,k u e* p K. (8) 

Here 7* is a virtual photon, and t = k* 2 . The branching ratio for the two-gamma decay 
mode is 

B(KL -> 11)= Mkiim t (9) 

v " ; 64tiT Kl ' v ; 



and the dilepton distribution of the Dalitz decays is given by 

fit) | 2 $(£, mi) T Kl ^, (10) 



dT KL ^ e +£~ J 



dt 

where f(t) = F(t)/F(0), me is the mass of the lepton, and 



a em 2 , i ,, . 2m» N / 4m? , , 
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There are two different sets of diagrams contributing to F(t). The first set are the pole 
terms shown in Figures 1 and 2 involving the virtual transitions Kl — > tt , ^s, ^1 followed 
by the two-photon decay modes driven by the WZ term. The second set of diagrams, shown 
in Figure 2, requires the weak interaction transitions Kl — ► 7r + 7r~7 and Kl — * 7r + 7r~77. 

As mentioned before, the 0(p A ) contributions due to the 7r° and 778 poles cancel each 
other due to the Gell-Mann-Okubo (GMO) relation ||. At 0(p 6 ) one-loop corrections 
and counterterms must be considered. From this various effects result. One is due to the 



deviation from the GMO relation, which is taken into account by replacing the lowest order 
masses by the physical masses. Another one is due to SU(3) breaking in virtual amplitudes 
Kl —> Vs with respect to Kl —* tt°, which we take into account by means of the parameter 
S below. A final effect stems from the one-loop corrections to the WZ term. As shown 
in flT5| , there is a t independent piece that is taken into account by simply replacing the 



decay constants in the chiral limit by the physical decay constants, while the t dependent 
piece requires renormalization. Fortunately, the latter only appears in contributions to the 
Kl — > £ + £~7 decay of order p 8 , and can therefore be altogether disregarded. The reason is 
that the ^-dependent corrections to the WZ term have the same relative weight for the 7r° 



and r] 8 as in the WZ term itself [15|1, and therefore the mechanism of cancellation by the 
GMO relation eliminates the order p 6 terms. Since the rji contribution is already of order 
p 6 at tree level, the one loop corrections are of order p 8 and, therefore, disregarded. From 
this, the pole piece of F(t) can be expressed as follows ||f|]: 

F! = rv + ^61 + ^82, (12) 

where 



r P = 1/ 1 
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®i = \ ((1 + 5) cosfl + 4/tsinfl) (-^ cos^-2 v / 2^ L sinfl 



VS M Vi 



B 2 = - ((1 + 5) sin6-4:Kcos6) ( ^- sin# + 2^ ^- cos9 ) . (13) 

We emphasize that F\ is very sensitive to the parameteres k and S, and therefore, cannot 
be reliably predicted. 

The one-loop diagrams of Figure 2 containing the vertices of order p 4 in (6) give the 
following contribution to F(t): 

^ 2 (t) = T |^{(a 2 + 2a 4 ) [lQ(Ml + Ml)-D(t^)]+CT(t^)} , (14) 



where CT denotes a counterterm, and the function D(t,n) is given in dimensional regular- 
ization by 



-10 ,, „ M 2 K , Ml 



D(t, /i) = £ [- + 2A - (log -f + log -£)] 

6 jJ, z jJ, z 

+A(F(M 2 n ,t) + F(Mlt)) , (15) 



where 



F( W ,^^(l- i )^log_ v _ H _ v _ 1 -2J m , 

y=-a, A = - + l + log47r-7 £; . (16) 

U-spin symmetry of £|A5=i| implies that the 7i~ + and the K + loops have the same sign. 

Only the ^-dependent piece needs renormalization provided by the counterterm CT(fi). 
The t-independent UV divergencies cancel when the two diagrams in Figure 2 are added. 
The t-dependent piece needs renormalization that is, as usual, accomplished by replacing in 
(14) CT(t,n) plus the term proportional to A by C(fi) t. As shown later, this counterterm 
can be estimated by analyzing the data. Therefore, ^(0) is unaffected by counterterm 
uncertainties, and is non- vanishing and proportional to M\ + M%. This term is thus a new 
addition to the Kl — > 77 amplitude. Notice that one can write a finite counterterm for 
the t = piece. Such a counterterm is included in F±, and in the pole model it is given 
by the 771 intermediate state as already discussed. The situation here resembles that of the 
strong interaction case at order p 4 , where the term proportional to the low energy constant 
L 7 only provides a finite renormalization, and is in fact dominated by the 771 pole ]I3]. The 



t independent pieces of F(t) vanish in the chiral limit in accordance with the Veltman- 
Sutherland theorem. 

The unknown parameters affecting the two-photon and Dalitz modes are JFi and (a 2 + 
204), while the counterterm C(fi)t only affects the Dalitz mode. When we neglect the 
t-dependence in F 2 , (10) gives 

ltU+^-7 = r ^-77 / 4m2 Ht, mi) dt (17) 



and the ratios with the experimental values are ||: 

-pExp 

i? e+e _ 7 = gp e+e ~ 7 = 0.97 ± 0.05 

-pExp 

R^i = l^ + ^ = 1.36 ± 0.14 (18) 

Thus, (17) is an excellent approximation in the e + e~ mode, where most of the rate is given 
by the low-t domain where the relative t-dependence is very small, while in the /i + /z~ case 
^min = 4m^ is large enough for the t-dependence to be noticeable. Thus, on the basis of the 
total rate of the /i + yU~ Dalitz mode it is clear that there is a significant t-dependence in F 2 , 



as was first experimentally noticed by the E799 Collaboration [11|. Denoting for the sake 
of convenience, 

C\ = (a 2 + 2a 4 ), 

c^+^Q-tof^fjj-m , , 19 , 

from the ratio -R^+^- 7 we find that they are approximately linearly related as 

C 2 ~ -sign(F(0)) (3.8 ± 1.3) + 1.1 C x . (20) 

With this relation we also find an excellent fit to both dilepton invariant mass 
distributions p|-[iT|. As expected, the Dalitz decays alone cannot pin-down C\ and the 



counterterm. It is necessary to invoke further observables for this to be possible. The 
value of C\ can also be estimated by means of two other processes, namely, Kl — ► 77 and 
Kl — > 7r + 7r~7. The amplitude of the first process is given in terms of F(0) = a 2 c ^p 8 L\ with 

L ^ = - F i + c ^ q^ F 2 ■ ( 21 ) 

From the Kl -^ 77 width (9) this combination has the value ±0.89. On the other hand, the 
second process permits a model dependent estimate of a different combination [12|], namely, 



L 2 = C X -F X . (22) 



While in (21) both terms are of the same chiral order, in (22) the first term is of lower 
order than the second. Indeed, C\ gives the strength of the direct emission Ml amplitude of 
order p A , while the second term appears in the correction of order p 6 to that amplitude ||12| . 



There is, in addition, the theoretically well known internal bremsstrahlung amplitude of 
order p 2 that is CP violating and cannot interfere with the direct emission one because it 
is of electric type. Thus, the direct emission decay rate can be cleanly identified, and the 
combination in (22) is then estimated to be in the interval 0.3 to 0.9 jn|. The positive 
sign of the combination is favored by a factorization model [12|, but the possibility of a 



negative sign is not altogether ruled out. For a similar negative interval the conclusions we 
draw here remain unchanged. Taking the estimated range at face value, we obtain C\ and 
F\ for the two possible signs of L\ as shown in the table. A first consequence is that the 
amplitude for the two-gamma decay has a potentially important addition beyond the pole 
terms, depending on the scenario, defined by the sign of L±, one considers. This is shown in 
the table by the ratio F 2 (0)/Fi, which can be as large as 0.6. 
The slope of the Dalitz decay form factor defined by 

'Ai^^^hrA^-^ x =w K < 23 > 

is consistent with the experimental value obtained from the e + e~ mode, &Exp = 0.6 ± 
0.25 H|l(|. In fact, in the acceptable range of C\ (20) gives b = 0.45 ± 0.2. We notice 



that the chiral logarithm terms as well as the counterterm give positive contributions to the 
slope. The fraction of the slope due to the counterterm is shown in the last entry of the 
table. We see that in one case it is dominating, giving support to the VMD model ||, while 
in the other case there are important chiral logarithm terms as well. 

Since the /x + /i~ mode shows more prominently the t-dependence, the results are mostly 
determined by that mode. In turn, for the e + e~ mode we predict i? e + e - 7 = 1.025 ± 0.010. 
This ratio is consistent within the 5% error of the experimental result (18). 

In conclusion, our leading order analysis shows two acceptable scenarios. In one of 
them there are sizeable long distance contributions from one-loop diagrams to both types of 



radiative decays considered, while in the other such contributions are small, implying that 
VMD models are an excellent picture. The question of which of the two scenarios is actually 
realized is clearly important. 
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FIGURES 
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Figure 1: Pole diagrams. The square represents the insertion 
of the | AS = 1 | effective non-leptonic weak interaction and the 
dot is the WZ term. 





Figure 2: Order p 6 one- loop diagrams. The square represents 
the insertion of £, 2 9 31 • The mesons in the loop are tt + and K + . 



TABLE 



u 


Ci 


Fi 


c 2 


^2(0)/^ 


b 


bcr/b 


0.89 
0.3 <L 2 < 0.9 


-0.6 ±0.6 


-1.2 ±0.3 


-4.5 ±1.3 


0.25 ±0.25 


0.5 ±0.2 


0.9 ±0.1 


-0.89 
0.3 <L 2 < 0.9 


3.1 ±0.6 


2.5 ±0.3 


7.2 ±1.3 


0.65 ± 0.15 


0.4 ±0.2 


0.4 ±0.1 



